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Abstract 

A search for charginos and neutralinos, predicted by supersymmetric theories, is per- 
formed using a data sample of 182.1 pb _1 taken at a centre-of-mass energy of 189 GeV with 
the OPAL detector at LEP. No evidence for chargino or neutralino production is found. 
Upper limits on chargino and neutralino pair production (x\Xi j X1X2) cross-sections are 
obtained as a function of the chargino mass ( m ^,± ) , the lightest neutralino mass ( m^o ) and 
the second lightest neutralino mass (m^o). Within the Constrained Minimal Supersym- 
metric Standard Model framework, and for m~± — m^o > 5 GeV, the 95% confidence level 
lower limits on m~± are 93.6 GeV for tan/3 = 1.5 and 94.1 GeV for tan/3 = 35. These 

Xi 

limits are obtained assuming a universal scalar mass mo > 500 GeV. The corresponding 
limits for all mo are 78.0 and 71.7 GeV. The 95% confidence level lower limits on the 
lightest neutralino mass, valid for any value of tan /3 are 32.8 GeV for tjiq > 500 GeV and 
31.6 GeV for all niQ. 
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1 Introduction 



Super symmetric (SUSY) extensions of the Standard Model predict the existence of charginos 
and neutralinos [|1]. Charginos, xf, are the mass eigenstates formed by the mixing of the fields 
of the fermionic partners of the charged gauge bosons (winos) and those of the charged Higgs 
bosons (charged higgsinos). Fermionic partners of the 7, the Z boson, and the neutral Higgs 
bosons mix to form the mass eigenstates called neutralinos, x® [] If charginos exist and are 
sufficiently light, they are pair-produced at LEP through 7- or Z-exchange in the s-channel. For 
the wino component, there is an additional production process through scalar electron- neutrino 
(z? e ) exchange in the t-channel. The production cross-section is large (several pb) unless the 
scalar neutrino (sneutrino) is light, in which case the cross-section is reduced by destructive 
interference between the s-channel and t-channel diagrams ||. In much of the parameter 
space, xt decays dominantly into Xi^ +hl or XiQcl' y i a a virtual W boson. For small scalar 
lepton masses, decays to leptons via a scalar lepton become important. R-parity conservation 
is assumed throughout this note. With this assumption, the Xi is stable and invisible Q and 
the experimental signature for X1X1 events is large missing momentum transverse to the beam 
axis. 

Neutralino pairs (X1X2) can produced through a virtual Z or 7 (s-channel) or by a scalar 
electron (t-channel) exchange 0]. The X2 wm decay into x*jW, Xi^ + ^~ or XiQQ through a 
virtual Z boson, sneutrino, slepton, squark or a neutral SUSY Higgs boson (h° or A ). The 
decay via a virtual Z is the dominant mode in most of the parameter space. For small scalar 
lepton masses, decays to a lepton pair via a scalar lepton are important. The experimental 
signature of X1X2 events is an acoplanar pair of leptons or jets. If the mass difference between 
X2 an d Xi is small, the experimental signature becomes monojet-like. 

Motivated by Grand Unification and to simplify the physics interpretation, the Constrained 
Minimal Supersymmetric Standard Model (Constrained MSSM) |], || is used to guide the 
analysis but the results are also applied to more general models. At the grand unified (GUT) 
mass scale, in the Constrained MSSM all the gauginos are assumed to have a common mass, 
m 1 /2i and all the sfermions have a common mass, m . Details are given in Section 5.2. 

Previous searches for charginos and neutralinos have been performed using data collected 
near the Z peak (LEP1), at centre-of-mass energies (y/s) of 130-136 GeV 0, 161 GeV 0, 
170-172 GeV @,|], and 183 GeV with luminosities of about 100, 5, 10, 10 and 60 pb" 1 



respectively. No evidence for signal has been found. 

In 1998 the LEP e + e~ collider at CERN was operated at \fs= 188.6 GeV. This paper reports 
on direct searches for charginos and neutralinos performed using the data sample collected at 
this centre-of-mass energy. The total integrated luminosity collected with the OPAL detector 
at this energy is 182.1 pb -1 . The selection criteria are similar to those used in ||10|| , but have 



been modified to improve the sensitivity of the analysis at the current energy. The description 



of the OPAL detector and its performance can be found in Ref. [12] and [ 13 



1 In each case, the index i = 1, 2 or j = 1 to 4 is ordered by increasing mass. 

2 The lightest supersymmetric particle (LSP) is either Xi or the scalar neutrino. We assume Xi is the LSP 
in the direct searches described in this paper. If the scalar neutrino is lighter than the chargino, xt ~> 
becomes the dominant decay mode. For this case, we use the results of Ref. || to calculate limits, as mentioned 
in Section 5.2. 
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2 Event Simulation 



Chargino and neutralino signal events are generated with the DFGT generator [14] which 



includes spin correlations and allows for a proper treatment of both the W boson and the 
Z boson width effects in the chargino and heavy neutralino decays. The generator includes 



initial-state radiation and uses the JETSET 7.4 package []IS| for the hadronisation of the quark- 
antiquark system in the hadronic decays of charginos and neutralinos. SUSYGEN p"6[ is used 



to calculate the branching fractions for the Constrained MSSM interpretation of the analysis. 

The sources of background to the chargino and neutralino signals are two-photon, lepton 
pair, multihadronic and four-fermion processes. Two-photon processes are the most important 
background for the case of a small mass difference between the Xi an d the Xi, or between the 
X® an d the x?) since such events have small visible energy and small transverse momentum. 
Using the Monte Carlo generators PHOJET [|T7|, PYTHIA fT§ and HERWIG JT§, hadronic 



events from two-photon processes are simulated in which the invariant mass of the photon- 
photon system is larger than 5.0 GeV. Monte Carlo samples for four-lepton events (e + e~e + e~, 

0- 



All other four- 



e + e~[x + fi~ and e + e - T + T~) are generated with the Vermaseren program 
fermion processes except for regions of phase-space covered by the two-photon simulations, 
are simulated using the grc4f generator f20fl , which takes into account all interferences. The 
dominant contributions come from W + W~, Wez/, 7*Z**) and ZZ^ processes. Lepton pairs 
are generated using the KORALZ [^lj] generator for t + t~(j) and /i + /i~(7) events, and the 
BHWIDE [22] program for e + e~ — > e + e~(7) events. Multihadronic (qq(7)) events are simulated 
using PYTHIA JTJ . 

Generated signal and background events are processed through the full simulation of the 
OPAL detector E3] and the same event analysis chain is applied to the simulated events as to 
the data. 



3 Analysis 

Calculations of experimental variables are performed as in ||. The following preselections are 
applied to reduce background due to two-photon events and interactions of beam particles with 
the beam pipe or residual gas: (1) the number of charged tracks is required to be at least 
two; (2) the observed transverse momentum of the whole event is required to be larger than 
1.8 GeV; (3) the energy deposited in each silicon-tungsten forward calorimeter and in each 
forward detector has to be less than 2 GeV (these detectors are located in the forward region, 
with polar angle f] | cos#| > 0.99, surrounding the beam pipe); (4) the visible invariant mass of 
the event has to exceed 2 GeV; (5) there should be no signal in the MIP plug scintillators 0. 



3 A right-handed coordinate system is adopted, where the x-axis points to the centre of the LEP ring, and 
positive z is along the electron beam direction. The angles 6 and <fi are the polar and azimuthal angles, 
respectively. 

4 The MIP plug scintillators [jl3| are an array of thin scintillating tiles with embedded wavelength shifting 
fibre readout which have been installed to improve the hermiticity of the detector. They cover the polar angular 
range between 43 and 200 mrad. 
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3.1 Detection of charginos 



The event sample is divided into three mutually exclusive categories, motivated by the topolo- 
gies expected for chargino events. Separate analyses are applied to the preselected events in 
each category: 

(A) iV c h > 4 and no isolated leptons, where iV c h is the number of charged tracks. The isolated 



lepton selection criteria are the same as those described in Ref. [10]. When both xt an d 



Xi decay hadronically, signal events tend to fall into this category for all but the smallest 
values of AM. (= m~+ — m^o). 

(B) iV c h > 4 and at least one isolated lepton. If just one of the xt decays leptonically, signal 
events tend to fall into this category. 

(C) iV c h < 4. Events tend to fall into this category if AM + is small or if both charginos decay 
leptonically. 

The fraction of xtxi events falling into category (A) is about 35-50% for most of the AM + 
range. This fraction drops to less than 15% if AM + is smaller than 5 GeV, since the average 
charged track multiplicity of the events becomes small. Similarly, the fraction of events falling 
into category (B) is also about 35-50% for most of the AM + range and is less than 10% if AM + 
is smaller than 5 GeV. In contrast, when AM + is smaller than 10 GeV, the fraction of events 
falling into category (C) is greater than about 50%, while if AM + is larger than 20 GeV, this 
fraction is about 10%. 

Since the chargino event topology mainly depends on the difference between the chargino 
mass and the lightest neutralino mass, different selection criteria are applied to four AM + 
regions: (I) AM+ < 10 GeV, (II) 10 GeV < AM+ < m-+/2, (III) m-+/2 < AM+ < m-+ - 
20 GeV, and (IV) m-+ —20 GeV< AM + < In region (I), background events come mainly 

from the two- photon processes. In regions (II) and (III), the main background processes are 
the four-fermion processes (W + W~, single W and 7*Z(*)). In these two regions the background 
level is modest. In region (IV) the W + W~ background becomes large and dominant. Since the 
W + W~ background is significant in the region of AM + > 85 GeV where the chargino decays 
via an on-mass-shell W-boson, a special analysis is applied to improve the sensitivity to the 
chargino signal for m~+ > 85 GeV and AM + > 85 GeV. Overlap between this analysis and the 
region (IV) standard analysis is avoided by selecting the analysis which minimises the expected 
cross-section limit calculated with only the expected number of background events. 

For each region a single set of cut values is determined which minimises the expected limit 
on the signal cross-section at 95% confidence level (C.L.) using the method of Ref. In this 
procedure, only the expected number of background events is taken into account and therefore 
the choice of cuts is independent of the number of candidates actually observed. 

The variables used in the selection criteria and their cut values are optimised in each AM + 



region. They are identical to the ones of Ref. [ID unless otherwise stated, and are therefore 
only briefly described in the following sections. 
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3.1.1 Analysis (A) (iV ch > 4 without isolated leptons) 



After the preselection, cuts on E{ w &/ E w - lS) |cos0 m i ss |, \P Z \ and \P z \/E v i S , are applied to further 
reduce background events from the two-photon and multihadronic processes. E vis is the total 
visible energy of the event, E iwd is the visible energy in the region of | cos^l > 0.9, 9 miss is the 
polar angle of the missing momentum and P z is the visible momentum along the beam axis. 
In regions (III) and (IV), the cut values for Ei wd /E w - ls and | cos^missl have been changed with 
respect to Ref. fllQ| , requiring them to be smaller than 0.15 and 0.90 respectively. Most of the 
two-photon background events are rejected by cuts against small p t HCAL and P t , the transverse 
momentum of the event measured with and without using the hadron calorimeter, respectively. 
In region (I), a cut is also applied against large P t HCAL (pHCAL <- gg Q e y) ^ Q rec l uce the W + W~ 
background. 



Jets are reconstructed using the Durham algorithm |25[ with jet resolution parameter y cut = 
0.005. By requiring the number of jets (iVj et ) to be between 3 and 5 inclusive, monojet events 
from the process "f*Z^ — > qquu are rejected for regions (I) and (II), and background events 
from qq(7) and We^ processes are reduced in regions (III) and (IV). 

In order to determine the acoplanarity angle (0 acO p), events passing the selection described 
above are forced into two jets, again using the Durham algorithm. Each jet is required to 
be far from the beam axis by cutting on its polar angle 6j (j = l,2). This cut ensures a good 
measurement of acop and further reduces the background from qq(7) and two-photon processes. 
The acoplanarity angle is required to be larger than 15° to reduce the qq(7) background. The 
acoplanarity angle distribution for region (III) is shown in Figure 1(a) before this cut. 

The cut on the visible mass, M vis , is optimised for each AM + region. If a lepton candidate 
(£') is found with an algorithm based on the "looser" isolation condition described in Ref. ||, 
the energy of this lepton, E#r, and the invariant mass calculated without this lepton, Mhad'j 
must be different from the values expected for W + W~ — > ivqq? events. The cuts on and 



M had / have been optimized with respect to Ref. |TUJ: for region (I) no cut is applied; for region 
(III) Mhad' is required to be smaller than 65 GeV and Eg smaller than 25 GeV. 

The background from W + W~ events and single-W events is efficiently suppressed by re- 
quiring that the highest and second highest jet energies, E\ and E 2 , be smaller than the typical 
jet energy expected for the W + W~ — > 4-jet process. The cuts on E\ and E 2 are identical to 
the ones in Ref. [|Kj] apart from region (II) where E\ is required to be between 2 and 30 GeV. 
In addition, in regions (III) and (IV), if E\ is larger than 40 GeV, M V - 1S is required to be ei- 
ther smaller than 70 GeV or larger than 95 GeV. In regions (III) and (IV) three-jet events 
with \P Z \ < 10 GeV are also rejected if M vis is close to the W mass. These cuts reduce the 
W + W~ — > ruqq' background with low-energy decay products of the r. 

A special analysis is applied in the region of AM + > 85 GeV, since the event topology of 
the signal is very similar to that of W + W _ — > 4 jets. After selecting well contained events 
with the cuts |cos# m i SS | < 0.95, E {wd /E vis < 0.15 and \P Z \ < 30 GeV, multi-jet events with 
large visible energy are selected with Aj et > 4 and 110 < E vis < 170 GeV. To select a clear 
4-jet topology 7/34 > 0.0075, 2/23 > 0.04 and 1/45 < 0.0015 are also required, where y/ n }{ n+ i} is 
defined as the minimum y cnt value at which the reconstruction of the event switches from n+1 
to n jets. Events having a "jet" consisting of a single 7 with energy greater than 20 GeV are 
considered to be 7qqg and are rejected. 

5 The acoplanarity angle, aC op, is denned as 180°- (j>, where (f> is the opening angle in the (r,4>) plane between 
the two jets. 
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The numbers of observed events and background events expected from the four different 
sources, for each AM + region and for the special analysis, are given in Table |I[ Typical detection 
efficiencies for XiXi events are 20-65% for AM + = 10-80 GeV in the standard analysis, and 
20-30% for m-+ > 85 GeV and AM+ > 85 GeV. 

Al 



Region 


I 


II 


III 


IV 


special 




AM + <10 GeV 


10 GeV<AM+ 


m_ + /2<AM+ 
x i 


m + -20 GeV 








<m_ + /2 


<m_+-20 GeV 


<AM+<m_+ 


AM+>85 GeV 


background 




77 


0.93 


0.46 


0.46 


0.46 


0.00 




0.00 


0.00 


0.02 


0.02 


0.00 


qq(7) 


0.07 


0.22 


0.83 


1.83 


4.57 


4f 


0.55 


1.27 


16.54 


24.66 


23.09 


total bkg 


1.5 ±0.7 


2.0 ±0.5 


17.8 ±0.9 


27.0 ± 1.1 


27.7 ± 1.0 


observed 


2 


3 


22 


29 


31 



Table 1: The numbers of expected background events for various Standard Model processes, 
normalised to the integrated luminosity, and the number of observed events in each 
AM_|_ region and for the special analysis for category (A). The errors in the total 
background include only the statistical error of the Monte Carlo samples. 

3.1.2 Analysis (B) (N ch > 4 with isolated leptons) 

After the preselection, cuts on | cos^missl and E^/ E v - ls are applied to reject two-photon back- 
ground events. The cut on E^/ E vis has been tightened with respect to Ref. fljj| : it is required 
to be smaller than 0.15 in regions (I) and (II) and smaller than 0.2 in regions (III) and (IV). 

In order to reject the W + W~ — > ivqq' background, the following cuts are applied: the 
momentum of isolated lepton candidates should be smaller than that expected from decays of 
the W (smaller than 15 GeV, 30 GeV, 40 GeV and between 5 and 40 GeV for regions (I) to 
(IV) respectively) and the invariant mass (M had ) of the event calculated excluding the highest 
momentum isolated lepton is required to be smaller than the W mass (smaller than 20 GeV and 
40 GeV for regions (I) and (II), between 10 and 60 GeV for region (III), and between 15 and 
70 GeV for region (IV)). The distribution of M had after the acO p cut is shown in Figure 1(b) 
for region (III). As is evident in this figure, most of the W + W~ background events are rejected 
by this cut. The invariant mass of the system formed by the missing momentum and the most 
energetic isolated lepton, M£ miss , is required to be larger than 110 GeV for region (IV). Finally, 
a M vis cut is applied to reject Wez/ events in which a fake lepton candidate is found in the 
W— > qq'(g) decay: it is required to be smaller than 30 GeV for region (I), between 25 and 
85 GeV for region (III), and no requirement on M vis is applied for region (IV). 

A special analysis is applied in the region of AM + > 85 GeV where there is a large W + W~ 
background. The selection criteria are identical to those in region (IV) up to the acO p cut. To 
reject further W + W~ — > iuqq' events while keeping a good signal efficiency, M had is required to 
be between 70 and 95 GeV, M vis between 90 and 125 GeV, and Me miss between 90 and 130 GeV. 

The numbers of observed events and background events expected from the four different 
sources, for each AM + region and for the special analysis, are given in Table [| Typical detection 



8 



efficiencies for XiXi events are 40-65% for AM + =10-80 GeV for the standard analysis, and 
20-30% for m-+ > 85 GeV and AM. > 85 GeV. 



Region 


I 


II 


III 


IV 


special 




AM+<10 GeV 


10 GeV<AM+ 


m_ + /2<AM + 


m.+ -20 GeV 

*l 








~ XT' 


<m +-20 GeV 

- *i 


<AM + <m + 


AM+>85 GeV 


background 




77 


3.80 


0.46 


0.00 


0.00 


0.00 


£+r( 7 ) 


0.02 


0.04 


0.11 


0.05 


0.02 


qq(7) 


0.00 


0.11 


0.14 


0.11 


0.04 


4f 


0.62 


1.38 


4.30 


6.37 


23.45 


total bkg 


4.4±1.3 


2.0±0.5 


4.6±0.4 


6.5±0.5 


23.5±0.9 


observed 


2 


1 


4 


7 


27 



Table 2: The numbers of expected background events for various Standard Model processes, 
normalised to the integrated luminosity, and the number of observed events in each 
AM+ region and for the special analysis for category (B). The errors in the total 
background include only the statistical error of the Monte Carlo samples. 



3.1.3 Analysis (C) (N ch < 4) 

This analysis is especially important for the region of AM + < 5 GeV. Because the background 
varies significantly with AM + in region (I), this region is split into 2 sub- regions (a,b). 

In order to reject events with charged particles which escape detection in the main detector, 
the net charge of the event is required to be zero. Since the signal is expected to have a two- 
lepton or two-jet topology, events are forced into two jets using the Durham jet algorithm |25|| . 
To improve the jet assignment, each jet must contain at least one charged track (A^ ch j > 1), 
must have significant energy (Ej >1.5 GeV) and the magnitude of the sum of the track charges 
(|Qj|) must not exceed 1. The Pt/E^ distributions for region (III) are shown in Figure 1(c) 
after these cuts. In region (I), if the acoplanarity angle is smaller than 70°, cuts are applied on 
P t , a t (the transverse momentum perpendicular to the event thrust axis), and | cos# a |, where 
9 a = tan _1 (a t /P z ). These cuts reduce the background contamination from two-photon and 
t + t~ processes. They are identical to those in Ref. |T(J for region (lb) but have been further 



optimized for region (la) where a t /E^ eSLm is required to lie between 0.025 and 0.1 and | cos 6^1 to 
be smaller than 0.9. To further reduce the two-photon background, cuts on P t and | cos^m^ 
are applied in region (I) if the acoplanarity angle is larger than 70° and in all other regions for 
any value of acoplanarity angle. The cuts on | cos^ m ; ss | are identical to those in Ref. [|IIJ] apart 
from region (la) where it is required to be smaller than 0.75. Pt/E hca _ m should be between 0.02 
and 0.04, between 0.03 and 0.05, and between 0.035 and 0.075 for regions (la), (lb) and (II), 
and should be larger than 0.095 and 0.1 for regions (III) and (IV). 

To reduce the background from e + e~/z + /i _ events in which one of the muons is emitted at 
a small polar angle and is not reconstructed as a good track, events are rejected if there is a 
track segment in the muon chamber or a hadron calorimeter cluster at a small polar angle, and 
within 1 radian in (r,0) of the missing momentum direction (P m i SS ). 
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Cuts on | cos 6^1 and aC o P are applied to reject two-photon, lepton-pair and 7*ZW — > l+t~vv 
events. The values of these cuts have been slightly modified with respect to Ref. []nj in that 
| cos#j| is now required to be smaller than 0.75 in region (la), and a cop should lie between 50° 
and 160° for region (la) and between 20° and 160° for region (lb). W + W~ — > l + vl~v events 
are rejected by upper cuts on M vis (20 GeV, 25 GeV, 30 GeV, 50 GeV and 75 GeV for regions 
(la) to (IV)) and on the higher energy of the two jets, Ei/E heam (identical to Ref. ]10 apart 
from region (la) where it is required to be smaller than 0.2). 

The numbers of background events expected from the four different sources, for each AM+ 
region, are given in Table |3|. The typical detection efficiencies for xfXi events are 30-65% for 
AM+ > 10 GeV, and the modest efficiency of 20% is obtained for AM+ = 5 GeV. 



Region 


I 


II 


III 


IV 


Sub-Region 


a 


b 


10 GeV<AM+ 


m_ + /2<AM + 
x i 


m +-20 GeV<AM+ 




3<AM + <5 GeV 


5<AAf+<10 GeV 


<m + /2 


<m + -20 GeV 


<m + 


background 




77 


6.89 


3.31 


9.07 


13.91 


12.99 




0.02 


0.16 


1.10 


2.97 


5.39 


qq(7) 


0.00 


0.00 


0.04 


0.04 


0.04 


4f 


0.95 


1.89 


9.72 


54.96 


97.97 


total bkg 


7.9±1.1 


5.4±0.8 


20.0±1.4 


71.9±2.1 


116.4±2.2 


observed 


12 


8 


23 


73 


112 



Table 3: The numbers of expected background events for various Standard Model processes, 
normalised to the integrated luminosity, and the number of observed events in each 
AM_|_ region for category (C). The errors in the total background include only the 
statistical error of the Monte Carlo samples. 



3.2 Detection of neutralinos 

To obtain optimal performance, the event sample is divided into two mutually exclusive cate- 
gories, motivated by the topologies expected for neutralino events. 

(C) iV c h < 4. Signal events in which decays into Xi^ + ^~ tend to fall into this category. 
Also, when the mass difference between x® and Xi (AM = m^o — m^o) is small, signal 
events tend to fall into this category independently of the x\ decay channel. 

(D) iV c h > 4. Signal events in which x% decays into XiQQ tend to fall into this category for 
modest and large values of AM . 

The event topology of X1X2 even ts depends mainly on the difference between the X2 and x± 
masses. Separate selection criteria are therefore used in four AM regions: (i) AM < 10 GeV, 
(ii) 10 < AM < 30 GeV, (iii) 30 < AM < 80 GeV, (iv) AM > 80 GeV. In regions (i) and (ii), 
the main sources of background are two-photon and 7*Z ( -*) — > qqz/P processes. In regions (iii) 
and (iv), the main sources of background are four-fermion processes (W + W~, Wez/ and 7*Z^*^). 
The fraction of events falling into category (C) is 10-20% for AM > 20 GeV but increases to 
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about 70% when AM < 5 GeV. The fraction of invisible events due to X2 ~^ Xi^ ~ y x\ vv 
decays is 20-30% depending on AM . 

The selection criteria applied for the low- multiplicity events (category (C)) in regions (i), 
(ii), (iii) and (iv) are identical to those used in analysis (C) of the chargino search in regions 
(la), (II), (III) and (IV), respectively (see Table [|). 

Events falling into category (D) have typically a monojet or a di-jet topology with large 
missing transverse momentum. The cuts described below are applied for these topologies. They 
are identical to the ones of Ref. pi unless otherwise stated. 



3.2.1 Analysis (D) (N ch > 4) 

To reduce the background from two-photon and qq(7) processes, cuts on | cos(6' miss )|, E iwd /E vis 
and missing transverse momenta are applied. The acoplanarity angle should be large to remove 
qq(7) background events. To ensure the reliability of the measurement of aC o P , both jets should 
have a polar angle 9j satisfying |cos#j| < 0.95. In region (iv), the aC op cut is loosened with 
respect to regions (i)-(iii), since the acoplanarity angle of signal events becomes smaller. To 
compensate for the resulting higher qq(7) background, the Ef wd /E v i S cut is tightened. 

After these cuts, the remaining background events come predominantly from 7*Z ( ** ) — > qquu, 
W + W~ — > ivqq and Wez/ — > qq eu. Cuts on the visible mass are applied to reduce W + W~ — > 
ivqq and Weu — > qq eu processes. 7*Z**) — > qquv background events are removed by a cut on 
the ratio of the visible mass to the visible energy, which is required to be larger than 0.3 for 
region (i) and 0.25 for regions (ii) and (iii). In regions (iii) and (iv), = 2/23-Evis < 30 GeV 2 is 
required to select a clear two-jet topology and to reject W + W~ — > ruqq' events. In Figure 1(d) 
the g?2 3 distribution is shown for region (iv) after all the other cuts. The numbers of background 
events expected from the four different sources, for each AMq region, are given in Table |4j. 
Typical detection efficiencies for X2X1 events are 50-65% for AM > 10 GeV. 



Region 


i 


ii 


iii 


iv 




AAf <10 GeV 


10 GeV<AAf <30 GeV 


30 GeV<AAf <80 GeV 


AAf >80 GeV 


background 




77 


3.06 


4.91 


1.48 


0.00 


£+£-(7) 


0.00 


0.07 


0.13 


0.33 


qq(7) 


0.00 


0.00 


0.00 


0.04 


4f 


0.25 


4.62 


11.72 


39.97 


total bkg 


3.3±1.1 


9.6±1.5 


13.3±1.0 


40.3±1.2 


observed 


2 


8 


14 


39 



Table 4: The numbers of expected background events for various Standard Model processes, 
normalised to the integrated luminosity, and the number of observed events in each 
AMq region for category (D). The errors in the total background include only the 
statistical error of the Monte Carlo samples. 
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4 Systematic uncertainties 



Systematic uncertainties on the number of expected signal and background events are estimated 
in the same manner as in the previous papers [§] and and are only briefly described here. 

For the number of expected signal events the uncertainties arise from the measurement of 
the integrated luminosity (0.5%), Monte-Carlo statistics and interpolation of the efficiencies 
to arbitrary values of m~± and m^o (2-10%), modelling of the cut variables in the Monte 
Carlo simulations (4-10%), fragmentation uncertainties in hadronic decays (< 2%) and detector 
calibration effects (< 1%). 

The angular distributions of the chargino and neutralino final-state decay products and their 
effect on the resulting signal detection efficiencies depend on the details of the parameters of the 
Constrained MSSM [![]. However, the corresponding variation of the efficiencies is determined 
to be less than 5% (relative), and this is taken into account in the systematic errors when 
obtaining the limits. Consequently, the limits are independent of the details of the Constrained 
MSSM. 

For the expected number of background events, the uncertainties are due to Monte Carlo 
statistics (see Tables 1 to 4), uncertainties in the amount of two-photon background (30%), 
uncertainties in the simulation of the four-fermion processes (17%), and modelling of the cut 
variables (< 7%), as determined in Ref. 0. 

In addition to effects included in the detector simulation, an efficiency loss of 2.9% (relative) 
arises from beam-related background in the silicon-tungsten forward calorimeter and in the 
forward detector which is estimated using random beam crossing events. 



5 Results 

No evidence for Xi~Xi and X1X2 production is observed. Exclusion regions and limits are 
determined by using the likelihood ratio method ||26|| , which assigns greater weight to the 
analysis which has the largest sensitivity. Systematic uncertainties on the efficiencies and on 
the number of expected background events are taken into account in the cross-section limit 



calculations according to Ref. [27 



5.1 Limits on the X1X1 an d X2X1 production cross-sections 

Figures 2(a) and (b) show model-independent upper-limits (95% C.L.) on the production cross- 
sections of X1X1 and X2X1, respectively. These are obtained assuming the specific decay 
mode X\ ~> xfW^*^ for XiXii and X2 ~^ Xi^*^ f° r X1X2 production. The results from 
the yfs =183 GeV [III analysis are also included in the limit calculation []. 



If the cross-section for X1X1 is larger than 0.75 pb and AM + is between 5 GeV and about 
80 GeV, m~+ is excluded at the 95% C.L. up to the kinematic limit. If the cross-section for 
X2X1 is larger than 0.95 pb and AM is greater than 7 GeV, m^o is excluded up to the kinematic 
limit at 95% C.L. 



6 In calculating limits, cross-sections at different y/s were estimated by weighting by (3/s, where (3 is 
Z±/Ebeam for xtxi production or p^o/E heam = p^o/E heam for X2X1 production. 
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5.2 Limits in the MSSM parameter space 



The phenomenology of the gaugino-higgsino sector of the MSSM is mostly determined by the 
following parameters: the SU(2) gaugino mass parameter at the weak scale (M2), the mixing 
parameter of the two Higgs doublet fields (//) and the ratio of the vacuum expectation values of 
the two Higgs doublets (tan/?). Assuming sfermions |§ and SUSY Higgs |28j sufficiently heavy 



not to intervene in the decay channels, these three parameters are sufficient to describe the 
chargino and neutralino sectors completely. Within the Constrained MSSM a large value of 
the common scalar mass, m (e.g., m = 500 GeV) leads to heavy sfermions and therefore to a 
negligible suppression of the cross-section due to interference with t-channel sneutrino exchange. 
Chargino decays would then proceed predominantly via a virtual or real W. On the other hand, 
a light tuq results in a low value of the mass of the sneutrino, enhancing the contribution of the 
t-channel exchange diagrams that have destructive interference with s-channel diagrams, thus 
reducing the cross-section for chargino pair production. Small values of mo also enhance the 
leptonic branching ratio of charginos. 

From the input parameters M 2 , fi, tan/5, m and A (the trilinear Higgs coupling), masses, 
production cross-sections and branching fractions are calculated according to the Constrained 
MSSM 0,^,0. For each set of input parameters, the total numbers of XiXii X2X1? X3X1 
and X2X2 even ts expected to be observed are calculated using the integrated luminosity, the 
cross-sections, the branching fractions, and the detection efficiencies (which depend upon the 
masses of the chargino, the lightest neutralino and next-to-lightest neutralino). Contributions 
from channels such as XiXa, X2X4; etc ■ ■ ■ are n °t included. The X3X1 channel is similar to the 
X2X1 channel, and cascade decays through X2 are taken into account. The relative importance 
of each of the analyses (A)-(D) changes with the leptonic or hadronic branching ratios, and 
the likelihood ratio method |26j is used to optimally weight each analysis depending on these 
branching ratios. 

Results are presented for two cases: (i) m = 500 GeV (i.e., heavy sfermions), and (ii) the 
value of mo that gives the smallest total numbers of expected chargino and neutralino events 
taking into account cross-sections, branching ratios, and detection efficiencies for each set of 
values of M 2 , fi and tan/?. This latter value of m leads to the most conservative limit at 
that point, so that the resulting limits are valid for all m . In searching for this value of m , 
only those values are considered that are compatible with the current limits on the sneutrino 
mass {rriy L > 43 GeV f29|), and upper limits on the cross-section for slepton pair production, 
particularly right-handed smuon and selectron pair production [[| . Particular attention is paid 
to the region of values of mo leading to the mass condition vrty « m-± by taking finer steps 
in the value of m . Note that we assume the stau mixing angle to be zero, and the added 
complication of possible enhanced decays to third-generation particles at large values of tan /3 
because of stau mixing is ignored. When m^ < m~±, resulting in a topology of acoplanar leptons 
and missing momentum, the upper limits on the cross-section for the two-body chargino decay 
from Ref. are used, while for m^ > m-± the three-body decays are dominant. Single photon 
topologies from X2X1 production and acoplanar photons with missing energy topologies from 
X2X2 w ith photonic decay X2 ~^ Xil are taken into account using the 95% C.L. cross-section 
upper limits on these topologies from OPAL results P0| . In both of these cases, if the relevant 



product of cross-section and branching ratio for a particular set of MSSM parameters is greater 
than the measured 95% C.L. upper limit presented in that paper, then that set of parameters 
is considered to be excluded. 
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The following regions of the Constrained MSSM parameters are scanned: < M 2 < 
2000 GeV, < 500 GeV, and A = ±M 2 , ±m and 0. The typical scan step is 0.2 GeV. 
Extensions beyond the scanned range have negligible effect on the quoted limits. The choice 
of A values is related to various scenarios of stop mixing, influencing the Higgs sector but 
having essentially no effect on the gaugino sector. No significant dependence on A is observed. 
Figure 3 shows the resulting exclusion regions in the (M 2 ,/i) plane for tan/5 = 1.5 and 35 with 
mo > 500 GeV and for all m . 

The restrictions on the Constrained MSSM parameter space presented here can be trans- 
formed into exclusion regions in the (m-±,m^o) or (m^o,m^o) plane. A given mass pair is 
excluded only if all considered Constrained MSSM parameters in the scan which lead to that 
same mass pair are excluded at the 95% C.L. The Xi mass limits are summarised in Table [5]. In 
the (m^±,m^o) plane, Figures 4(a) and (b) show the corresponding 95% C.L. exclusion regions 
for tan/5 = 1.5 and 35. Figures 4(c) and (d) show the corresponding 95% C.L. exclusion regions 
in the (m^o,m^o) plane, for tan/5 = 1.5 and 35. Mass limits on Xii X21 an d Xs are summarised 
in Table []. 







tan/5 =1.5 


tan /5 = 35 


m > 500 GeV 


AM+ > 5 GeV 


m-+ > 93.6 GeV 


m-+ > 94.1 GeV 


All mo (see text) 


AM+ > 5 GeV 


m*+ > 78.0 GeV 


m.+ > 71.7 GeV 



Table 5: Lower limits at 95% C.L. obtained on the lightest chargino mass. 







tan/5 = 1.5 


tan/5 = 35 


m > 500 GeV 


No AMo restriction 


m^o > 40.2 GeV 


m^o > 48.5 GeV 


AM > 10 GeV 


m^. > 67.8 GeV 
m^> 106.0 GeV 


m,o > 94.3 GeV 

A2 

m^o > 124.0 GeV 


All mo (see text) 


No AMo restriction 


m^o > 34.1 GeV 

Al 


m^o > 38.9 GeV 


AM > 10 GeV 


m^o > 55.9 GeV 
mJ> 106.6 GeV 


m^o > 86.2 GeV 

A2 

m^o > 134.3 GeV 

A3 



Table 6: Lower limits at 95% C.L. obtained on m^-o, m^o, and m^o. 



Figure 5 shows the dependence of the mass limits on the value of tan/5. Of particular 
interest is the absolute lower limit, in the framework of the Constrained MSSM, on the mass 
of the lightest neutralino of m^o > 32.8 GeV (31.6 GeV) at 95% C.L. for m > 500 GeV (all 
mo). This has implications for direct searches for the lightest neutralino as a candidate for dark 
matter [01 ]. Since the formulae for couplings and masses in the gaugino sector are symmetric 



in tan/5 and 1/ tan/3, these results also hold for tan/5 < 1. 



6 Summary and Conclusion 

A data sample corresponding to an integrated luminosity of 182.1 pb -1 at ^/s =188.6 GeV, 
collected with the OPAL detector, has been analysed to search for pair-production of charginos 
{XiXi) an d neutralinos (X2X1) predicted by supersymmetric theories. Decays of Xi Xi^ ±iy 
or Xi°Sl' an d decays of x 2 mto x'jW; Xi^ + ^~ or XiQQ are looked for. No evidence for such events 



14 



has been observed. The exclusion limits on xf an d X° production are significantly higher with 
respect to the results obtained at lower centre-of-mass energies. 

Exclusion regions valid at the 95% confidence level have been derived in the framework of 
the Constrained MSSM, in which only three parameters (M2, fi and tan/3) are necessary to 
describe the chargino and neutralino sectors. These restrictions in parameter space have been 
transformed into mass limits valid at the 95% confidence level. Assuming m ~± — m^o > 5 GeV, 
the lower mass limit of the chargino is 93.6 GeV for tan/3 = 1.5 and 94.1 GeV for tan/3 = 35 
for the case of a large universal scalar mass (m > 500 GeV); for all m , the mass limit is 
78.0 GeV for tan/3 = 1.5 and 71.7 GeV for tan/3 = 35. The lower mass limit for the lightest 
neutralino is 32.8 GeV for the case of itiq > 500 GeV and 31.6 GeV for all m . This latter 
result has implications for searches for the lightest neutralino as a dark matter candidate. 
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Figure 1: Distributions of some essential observables used to select chargino and neu- 
tralino events: (a) acoplanarity angle in analysis (A), (b) invariant mass of the event 
calculated excluding the highest momentum isolated lepton in analysis (B), (c) ratio of 
transverse momentum to beam energy, P t /E\>, in analysis (C), and (d) d| 3 in analysis 
(D). The data are shown with error bars and the distributions from the background 
processes are shown as filled histograms: dilepton events (double hatched area), two- 
photon processes (negative slope hatching area), four-fermion processes including W-pair 
events (positive slope hatching area), and multihadronic events (open area). The arrows, 
pointing into the accepted regions, show where the analysis cuts are applied. In Figures 
(a)-(c) the dashed line shows the prediction for a chargino signal with = 94 GeV and 
m^o = 47 GeV. In Figure (d) the dashed line shows the prediction for a neutralino signal 
with m^o = 145 GeV and m^o = 35 GeV. The normalisations of the signal histograms are 
arbitrary. 
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Figure 2: Contour plots of the 95% CL. upper- limits on the production cross-sections 
at ,/s = 189 GeV (a) cr~+~- for e + e~ — > xtx 



Xi is assumed to decay only into XiW^ , and the 



XiXl- The 



and (b) cr^o^o for e + e 

is assumed to decay only into 
Xi^*\ The cross-hatched region, for which m^o + m^o < mz, is covered by searches 
at LEP1, which set upper-limits on cr^o^o of order 5 pb, and is not considered in this 
analysis. The kinematical boundaries for XiXi an d X1X2 are shown as dashed lines. If 

there is no exclusion region in the (m-+,m^o) plane. If 

85 GeV 
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in the (m^o,m^o) plane. If c^o^o is between 0.11 and 0.13 pb a small region is excluded 



CJ-+-- is smaller than 0.068 pb, tnere is no exclusion region m tne [m-+,m 9 o 
CT-+-- is between 0.068 and 0.075 pb a small region is excluded around m-4 
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Figure 3: Exclusion regions at 95% C.L. in the \x-Mi plane within the framework of 
the Constrained MSSM for (a) tan/3 = 1.5 and (b) tan/3 = 35. The speckled areas 
show the regions excluded by LEP1 data, and the shaded areas the additional exclusion 
regions using the data from ^Js = 183 and 189 GeV. The dark-shaded regions are valid for 
all values of m , and the light-shaded regions are the additional excluded parameters if 
mo > 500 GeV (i.e., heavy scalar leptons). The kinematical boundary at ^fs = 189 GeV 
for XiXi production is shown by dashed lines. 
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Figure 4: The 95% CL. excluded regions in the (m^+ ,m^o) plane within the framework of the 
Constrained MSSM for the case of any value of (dark-shaded regions) and itlq > 500 GeV 
(extending to light-shaded regions) for (a) tan/3 = 1.5 and (b) tan/3 = 35. The thick solid lines 
represent the theoretical bounds of the Constrained MSSM parameter space. The kinematical 
boundaries for production and decay at \/s = 189 GeV are shown by dashed lines. Similar 
plots for neutralino masses are shown in (c) for tan/3 = 1.5 and (d) for tan/3 = 35. 
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Figure 5: The 95% CL. mass limit on the lightest neutralino Xi as a function of tan/3 within 
the framework of the Constrained MSSM. Exclusion regions for all tuq values (dark-shaded 
region) and m > 500 GeV (extending to the light-shaded region). 
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